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ABSTRACT 

This is the second paper of a series reporting the results from the PopStar evolutionary 
synthesis models. Here we present synthetic emission line spectra of Hii regions photoionized 
by young star clusters, for seven values of cluster masses and for ages between 0.1 and 5.2 
Myr The ionizing Spectral Energy Di s tributi ons (SEDs) are those obtained by the PopStar 
code dMolla. Garcia- Vargas & BressanL l2009h for six different metallicities, with a very low 
metallicity set, Z=0.0001, not included in previous similar works. We assume that the radius 
of the Hii region is the distance at which the ionized gas is deposited by the action of the 
mechanical energy of the winds and supernovae from the central ionizing young cluster In 
this way the ionization parameter is eliminated as free argument, since now its value is ob- 
tained from the cluster physical properties (mass, age and metallicity) and from the gaseous 
medium characteristics (density and abundances). We discuss our results and compare them 
with those from previous models and also with a large and data set of giant Hii regions for 
which abundances have been derived in a homogeneous manner The values of the [OIII] lines 
(at AA 4363, 4959, 5007 A) in the lowest metallicity nebulae are found to be very weak and 
similar to those coming from very high metallicity regions (solar or over-solar). Thus, the sole 
use of the oxygen lines is not enough to distinguish between very low and very high metal- 
licity regions. In these cases we emphasize the need of the additional support of alternative 
metallicity tracers, like the [SIII] Unes in the near-IR. 

Key words: ISM: HII regions - ISM: abundances - ISM: evolution - ISM: hnes and bands 



1 INTRODUCTION 

Hn regions have been wi del y stud ied for the last thre e d ecades 



I Pagel. Edmunds & SmitH, 1 198C 



Dopita & Evans . 1986; Diaz et al. 



Evans & Dopital 



1 19871 : IVflchez et al.L 



1985 



198S 



Diaz et al.l Il99ll ; Ipfai Il994l) . In those studies the functional 
parameters, such as ionization parameter, effective temperature 
of the ionizing stars and oxygen abundance, were translated into 
physical parameters associated to the star clusters (mass, age and 
metallicity). In order to do this, several grids using evolutionary 
models plus photoioniz at ion co d es were computed by our group 
dOarcfa Vargas & Dfazl 1 19941: iGarcia- Vargas. Bressan & Dfazl 



Il995a|| h 



Garcia- Vargas. Molla & Bressanl 119981) and oth- 



Stasinskal 1 19801 : IStasinska et al.L Il98ll : 



ers dStasinskal 1978 
IStasinskalll98l ll99C 

In the last decade the evolutionary synthesis technique has 
been improved with the inclusion of updated stellar tracks or 
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isochrones, the use of better stellar spectra as input for the codes, 
a better treatment of binary stars, a wider wavelength coverage 
and a higher spectral resolu t ion. C odes like starburst99 (here- 
inafter stb99 iLeitherer et al.L 1 19991) , use the Geneva group stel- 
lar tracks and are t uned for star forming regions. On the contrary 
codes like pegase jpioc & Rocca-Volmerang3 . 1 19971) or galaxev 
teruzual & CharlotlllOoi) ^ are specially computed for intermediate 
and evolved populations and use the Padova group stellar tracks. 
These codes, among others, have been intensively applied to spec- 
trophotometric catalogues to derive the physical properties of stel- 
lar populations. 

The results of these codes, mainly the SEDs, have been 
used in turn by many authors to compute models for star form- 
ing regions, Hn and starburt galaxies, or have been applied 
to particular cases. Recently, a more realistic approach to the 
Hn region structure, and therefore the ionization bubble, has 
been considered to improve pho t oioniz ation model results. Thus, 
IMov. Rocca-Volmerange & Fiod ( 1200 ih analyzed coherently the 
stellar and nebular energy distributions of starbursts and Hn galax- 
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ies, by using p e gase i n addition to the photoionization code cloudy 
jperland et all 1 19981) . constructing models in which the filling 
factor and the radius of the nebula are fixed to obtain a range 
of values of the ionization paramete r similar to those derived 
from observations. iDopita et al.l j2000l) used both pegase and stb99 
models to compute the SEDs of young star clusters. Main dis- 
crepances between the SEDs of these two grids of models were 
driven by the dilferent evolutionary tracks used, Padova's for pe- 
gase and Geneva's for stb99. They studied the emission line se- 
quence and computed the Hii region spectra as a function of age, 



value the young cluster mass using dust-free models. The dust re- 
emits the light in the mid and far infrared spectral regions, with 
evident conseque nces for the population synthesis in these spec- 
tral windows iBressan, Granato & Silvall 19981: ISilva et all 199^ 



metallicity and i oniza tion parameter usi ng the photoionization code 
MAPPINGS vO.3 ([Suthe rland & Dopitj 1 19931 : iDopita etafl |2002|: 
lOroves. Dopita & Sutherland! 120041) . More recently lOopita et al.l 



l 2006ij) used stb99 plus mappings III to produce a self-consistent 
model. In this model they consider that the expansion and internal 
pressure of the Hn regions depends on the mechanical energy from 
the central cluster, replacing the ionization parameter by a more re- 
alistic one dependent on cluster ma ss and pressure in the ISM . A 
similar approach has been used by IStasinska & Izotovl ( |2003|) . In 
this work they calculate photoionization models of evolving star- 
bursts by using as ionizi ng continuum the SEDs from the code by 
ISchaerer & Vaccal ( 1 199^ ). These models were computed for differ- 
ent metallicity bins to compare with the con'esponding metallic- 
ity observations, assuming the geometry of an adiabatic expand- 
ing bubble, instead of that of a region with a fixed radius (like in 
iGarcia- Vargas et ali[T9953lbl : IStasihska & Leithere3.ll996h . 

The present work is the second paper of a series of 
three dedicated to the PopStar models description and initial 
test-cases application. PopStar is a new grid of evolutionary 
synthe sis models described in iMolla. Garcia- Vargas & BressanI 
( l2009l ,hereinafter Paper I), where their suitability to model 
stellar populations in a wide range of ages and metal- 
lic ities is shown. T h ese m o dels are an updated version 
of Garcia- Vargas et al. 1 1995bV 'Garcia- Varg as. Molla & BressanI 
( 1 19981) and iMolla & Gar cia- Vargas ( 2000). One of the main 
improvements of this new grid is a very careful treatment 
of the emerging spectrua of hot stars, which can be ei- 
ther massive stars or post-AGB objects. In previous works 
by our group dGarcfa Vargas & Dfa 3, 1 1994 [g arcia- Vargas et all 
Il995al lbl) h ot massive stars we re mod elled using the atmo- 
spheres of IClegg & Middlemasj h987h . In the new PopStar 
SEDs, NLTE blanketed models for 0,B and WR stars by 
ISmith. Norris & CrowtheJ j2002h have been used. These models 
include a detailed treatment of stellar winds which modifies the fi- 
nal SEDs obtained for young stellar populations and consequently 
the emission line spectrum of the surrounding nebula. 

Our aim in the present work is to compute a set of updated 
photoionization models, and therefore emission line spectra, for Hii 
regions following the evolution of the ionizing young star cluster, 
whose SED is computed by the evolutionary synthesis models Pop- 
Star. We also compare these results with a sample of HII regions 
where the metallicity has been carefully derived through the use of 
appropriate calibrators. 

This grid of models has been computed without taking into 
account the presence of dust and this imposes some limitations to 
its use. Our purpose is to compute a grid of models covering the 
range of physical parameters found in GEHR where the dust effect 
is weak on the optical spectral lines, which are the ones computed 
in this work. However, in starburt galaxies, where star formation 
is taking place inside molecular clouds where the dust is present, 
its effect can be very important. The dust will absorb a fraction of 
the emitted ionizing photons and therefore will affect the derived 



IPanuzzo et"all |2003| ; ISressan et all l2006l: IClemens et all (2009t) 
and even in the radio wavelength range dVega et alll2008l) . Our pur- 
pose is to compute in the near future self-consistent models includ- 
ing the chemical and the spectro-photometric evolution, for spiral 
and irregular galaxies, where star formation and dust effects are im- 
portant, (since star formation takes place inside molecular clouds). 
In this case, the molecular clouds will be specifically treated as a 
separated phase, the enrichment of the gas will be well followed 
and the dust will be included consistently. 

Finally, we must mention that one of the main differences 
between the isochrones from Geneva and Padova groups is the 
inclusion of stellar rotation. Stellar rotation affects the evolu- 
tionary tracks, lifetimes, and chemical compositions of massive 
stars, as well as the formation of red super-giants and WR stars 
dMaeder & MevneA[20o3 ; lMeynet & MaedeJ.l2005l) . In particular, 
the rotation of a star plays a determinant role at very low metallic- 
ities, producing high mass loss where almost none was expected, 
as explained by Ekstrom et al. (2008, and references therein). The 
revised grid of stellar evolutionary tracks accounting for rotation, 
recently released by the Geneva group, has been implemented into 
the stb99 evolutionary s ynthesis code in the preliminary models 
bv lVazQuez etlD JioOTh. Massive stars are predicted to be hot- 
ter and more luminous than previously thought, this effect being 
higher for decreasing metallicity. Individual stars now tend to be 
bluer and more luminous, increasing by a factor of 2 (or even more) 
the light-to-mass ratios at ultraviolet to near-infrared wavelengths, 
as well as the total number of ionizing photons. However, we have 
not taken this effect into account because we have used SEDs based 
in a set of revised isochrones from the Padova group. This revision 
affects mainly the computation of intermediate age stellar popula- 
tions (RGB and ABB phases), and does not include rotation. 

Section 2 gives a summary of our theoretical evolutionary syn- 
thesis and photoionization models, the inputs and the hypotheses 
assumed for the grid calculation. Section 3 shows our model re- 
sults in terms of the time evolution of emission lines. In this sec- 
tion we also include the co mparison with the previous models by 
iGarcia- Vargas et alj ( ll995a llbl). In Section 4 we present a discussion 
of the models based on the comparison with observed emission line 
ratios which have been carefully compiled from the literature for 
objects whose abundances have been consistently re-calculated us- 
ing the appropriate techniques. This section also includes the guide- 
lines for deriving the physical properties of the ionizing star clusters 
through the analysis of the correponding gas emission line spectra. 
Finally, our conclusions are given in Section 5. 



2 SUMMARY OF THEORETICAL MODELS 
2.1 Evolutionary synthesis models. 

The SEDs used in this work have been taken from the new Pop- 
Star evolutionary models (Paper I). The basic model grid is com- 
posed by Single Stellar Populations (hereinafter SSPs) for 6 dif- 
ferent Initial Mass Functions (IMF) from which only those with a 
Salpeter power law {Salpeter, \95% and different mass limits have 
been used here: 0.85Mo and 120Mo (hereinafter SALl), O.ISMq 
and IOOMq (hereinafter SAL2), and IMq and IOOMq ( hereinafter 
STB). None of the models include either binaries or mass segrega- 
tion. 
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Isochrones are an update of those from iBressan et al" I l ll998l) 
for 6 different metallicities: 0.0001. 0.0004, 0.004, 0.008, 0.02 and 
0.05 the lowest metallicity set being now added. The age cover- 
age is from log/ = 5.00 to 10.30 with a variable time resolution 
of A(logf) = 0.01 in the youngest stellar ages. The details of the 

isochrones are described in Paper I. 

Atmosphere models are from iLeieune. Cuisinier & BuseJ 

( 1 19971) with an excellent coverage in effective temperature, grav- 
ity and metallicities, for stars with Teff < 25000K. For O, 
B and WR the code uses th e NLTE blanketed models by 
ISmith. Norris & CrowtheJ | |2002|) at Z = 0.001, 0.004, 0.008, 0.02 
and 0.04. There are no models available for Z = 0.0001, 0.0004 
and Z=0.04. For these stars, we have selected the models with 
the closest metallicity, that is, Z=0.04 for the isochrones with 
Z=0.05 and Z=0.001 for isochrones with Z=0.0001 and 0.0004. 
The possible misselection for the lowest metallicities is not im- 
portant for what refers to WR stars since, in principle, there are 
no WR at these lowest Z. This is in any case a limitation of the 
models: the grid for NLTE atmosphere models is not as fine as 
the LTE Lejeune. There are 110 models fo r 0-B stars, calculated 
by IPauldrach. Hoffmann & LennorJ ( 1200 ih . with 25000K < Tefi' 

< 51500K and 2.95 < log;; < 4.00, and 120 mode ls for WR stars 
(60 WN and 60 WC), from lffillier & Millej l ll998l) . with 30000K 

< r* < 120000K and 1.3i?o < R' < 20.3i?o for WN, and with 
40000is: < r* < 140000K and O.SRo < R' < 93Ro for WC. T* and 
R* are the temperature and the radius at a Roseland optical depth 
of 10. The assignation of the appropriate WR model is consistently 
made by using the relationships among opacity, mass loss and ve- 
locity wind described in Paper I. 

For post-AGB and PN wit h Tefi^ higher t han 50000K and up to 
220000K the NLTE models by!Rauch| l l2003h are taken. For higher 
temperatures PopStar uses black bodies. The use of these last mod- 
els affects the resulting intermediate age SEDs, which are not used 
in the present work. 

2.2 Photoionization models for Hn regions. 

We have studied the evolution of a cluster along the first 5.2 Myr, in 
21 steps of time. Clusters older than this age do not have the neces- 
sary ionization photons to produce a visible emission line spectrum, 
although they are still detectable on H„ imagefl Ionizing clusters 
have been assumed to form in a single burst with different masses. 
We have run models assuming seven values for the total cluster 
mass, 0.12, 0.20, 0.40, 0.60, 1.00, 1.50 and 2.00 xlO' Mq, using 
SAL2 IMF (with m;o„,=0.15Mo and m„p=100 Mq). Masses in this 
range are able to provide the observed number of ionizing photons 
of most medium to large extragalactic Hii regions. Also, in order to 
compare with our old photoionization models we have used SALl 
IMR 

We hav e use d the photoionization code CLOUDY 
jFerland" et al.L 1 19981) to obtain the emission line spectra of 
the modelled Hn regions for different metallicities. The gas is 
assumed to be ionized by the massive stars of the young cluster 
whose ionizing spectra have been taken fro the PopStar SEDs 



See Garcia- Vargas, Moffa & Martm-Manjon, in preparation, where we 
calculate the photometric properties of young star clusters taking into ac- 
count the contamination by the emission lines and the underiying old pop- 
ulations where these clusters are embedded in the computed colors, for an 
analysis about the cluster physical properties determined from the photo- 
metric information only 



Table 1. Abundances, logX/H used in the models jGrevesse & Sauval 
1 19981) 



El. 


sol. 


sol. 


0.0001 


0.0004 


0.004 


0.008 


0.02 


0.05 






depl 


depl 


depl 


depl 


depl 


depl 


depl 


He 


-1.09 


-1.09 


-1.24 


-1.24 


-1.20 


-1.16 


-1.07 


-0.90 


C 


-3.48 


-3.48 


-5.71 


-5.11 


-4.11 


-3.81 


-3.41 


-3.01 


N 


-4.08 


-3.96 


-6.19 


-5.59 


-4.59 


-4.29 


-3.89 


-3.49 


O 


-3.17 


-3.18 


-5.41 


-4.81 


-3.81 


-3.51 


-3.11 


-2.71 


Ne 


-3.92 


-3.92 


-6.15 


-5.55 


-4.55 


-4.25 


-3.85 


-3.45 


Na 


-5.67 


-6.67 


-8.90 


-8.30 


-7.30 


-7.00 


-6.60 


-6.20 


Mg 


-4.42 


-5.42 


-7.65 


-7.05 


-6.05 


-5.75 


-5.35 


-4.95 


Al 


-5.53 


-6.53 


-8.76 


-8.16 


-7.16 


-6.86 


-6.46 


-6.06 


Si 


-4.45 


-4.75 


-6.98 


-6.38 


-5.38 


-5.08 


-4.68 


-4.28 


S 


-4.67 


-4.67 


-6.90 


-6.30 


-5.30 


-5.00 


-4.60 


-4.20 


Ar 


-5.60 


-5.60 


-7.83 


-7.23 


-6.23 


-5.93 


-5.53 


-5.13 


Ca 


-5.64 


-6.64 


-8.87 


-8.27 


-7.27 


-6.97 


-6.57 


-6.17 


Fe 


-4.50 


-5.50 


-7.73 


-7.13 


-6.13 


-5.83 


-5.43 


-5.03 


Ni 


-5.75 


-6.75 


-8.98 


-8.38 


-7.38 


-7.08 


-6.68 


-6.28 



computed in Paper I, as explained in the previous section. The 
SEDs are given for a normalized mass of 1 Mq, hence the stellar 
mass of the cluster must be used to scale the number of ionizing 
photons and other absolute parameters given by PopStar for each 
IMF. The shape of the ionizing continuum, the number of ionizing 
photons, Q(H), and the Hn region radius, R^, set by the action 
of the cluster mechanical energy, are obtained directly from the 
ionizing cluster parameters: mass, age and metallicity. 

The chemical composition of the gas and its spatial distribu- 
tion, assuming a given geometry, together with the medium den- 
sity, are used as inputs for the photoionization code. It is as- 
sumed that both the cluster and the surrounding gas have the 
same chemical compositio n. The solar abundances are taken from 
iGrevesse & Sauval ( Il998h . The use of these solar abundances im- 
plies that Zq = 0.017 Q thus, Z = 0.02 of our models does not 
correspond to the solar value but to 1.17 Zq. 

The solar abundances are summarized in Table[T] Column 2. 
Some refractary elements, Na, Al, Si, Ca, Fe and Ni, must be de- 
pleted due to the presence of dust grains mixed with the ionized 
gas which can affect the UV photon absorption and decrease the 
electronic temperatu re. Solar depleted abundances, obtained using 
lOamettetaf] ( 119951) depletion factors, are in Column 3. For each 
metallicity, the element abundances heavier than helitmi have been 
scaled by a constant factor, with respect to the hydrogen content, 
according to the solar depleted abundances, as given in Columns 4 
to 8 of the table. 

The radius of the modelled region has been derived from 
the mechanical energy produce d by the expanding atmospheres of 
massive stars with strong winds. ICastor. McCrav & Weave3 ^19751) 
demonstrated that an early-type star with a strong stellar wind can 
blow out a large cavity or bubble in the surrounding gas. The wind- 
driven shell begins to evolve with an initial phase of free expansion 
followed by an adiabatic expansion phase, and then the material 
shall collapse into a thin, cold shell as a result of radiative cool- 
ing. At this stage the gas shall trap the ionization front and thus 
the radiative phase begins. In this phase the ionizing photons are 



^ Asplund, Grevesse & Sauval (2005) obtained a solar abundance Zq = 
0.012. This value is still qu estioned I'B ahcall, Serenelli & Basul |2005| : 
brake & Test3 . l2005l : lAntia & Basu , .2005,) because it does not fit the he- 
lioseismological constraints. 
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absorbed and the region cools via emission in tlie Balmer lines. In 
this process, the radius of the outer shock, R,, evolves as: 

= l.6(6/n)"'P"(pc) (1) 
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Table 2. Emission line spectrum from a Hn region of mass 4x10" Mo and Z=0.008 as a function of cluster parameters for n;/ = 10 cm . This is Table 2a. 
Table 2b includes the models with n/y = 100 cm"'. Both tables for the whole set of cluster masses and metallicities are provided in electronic format. 

Zmet logage Mass [Oil] [QUI] [GUI] [GUI] [OT] [Nil] [Nil] [SlT] 
yr Mo 3727A 5007A 4959A 4363A 6300A 6584A 6548A 6716A 



0.008 


5.00 


410" 


0.7398 


3.5490 


1.1791 


0.0140 


0.0189 


0.2101 


0.0712 


0.1413 


0.008 


5.48 


410" 


0.7644 


3.8260 


1.2711 


0.0162 


0.0206 


0.2121 


0.0719 


0.1492 


0.008 


5.70 


410" 


0.7802 


3.9562 


1.3144 


0.0173 


0.0214 


0.2147 


0.0728 


0.1535 


0.008 


5.85 


410" 


0.8049 


4.1260 


1.3708 


0.0186 


0.0222 


0.2189 


0.0742 


0.1585 


0.008 


6.00 


410" 


0.8196 


3.8744 


1.2872 


0.0168 


0.0221 


0.2246 


0.0761 


0.1592 


0.008 


6.10 


410" 


0.8279 


3.6580 


1.2153 


0.0154 


0.0217 


0.2282 


0.0773 


0.1583 


0.008 


6.18 


410" 


0.8788 


3.5548 


1.1810 


0.0147 


0.0221 


0.2386 


0.0809 


0.1638 


0.008 


6.24 


410" 


1.0206 


3.3794 


1.1227 


0.0135 


0.0211 


0.2476 


0.0839 


0.1671 


0.008 


6.30 


410" 


1.2825 


2.9609 


0.9837 


0.0111 


0.0201 


0.2722 


0.0923 


0.1745 


0.008 


6.35 


410" 


1.4038 


2.4940 


0.8286 


0.0086 


0.0198 


0.3065 


0.1039 


0.1884 


0.008 


6.40 


410" 


1.6003 


1.8690 


0.6209 


0.0057 


0.0182 


0.3542 


0.1200 


0.2073 


0.008 


6.44 


410" 


1.6950 


1.6437 


0.5461 


0.0048 


0.0183 


0.3994 


0.1354 


0.2296 


0.008 


6.48 


410" 


1.7781 


1.2606 


0.4188 


0.0032 


0.0164 


0.4705 


0.1594 


0.2616 


0.008 


6.51 


410" 


1.4910 


1.5379 


0.5109 


0.0042 


0.0221 


0.4548 


0.1541 


0.2791 


0.008 


6.54 


410" 


1.8240 


2.2962 


0.7629 


0.0082 


0.0459 


0.5305 


0.1798 


0.3893 


0.008 


6.57 


410" 


2.2596 


2.5285 


0.8400 


0.0099 


0.0712 


0.6845 


0.2320 


0.5287 


0.008 


6.60 


410" 


2.8121 


2.1852 


0.7260 


0.0083 


0.1159 


0.9505 


0.3221 


0.7827 


0.008 


6.63 


410" 


2.2239 


1.3071 


0.4342 


0.0037 


0.0838 


0.8983 


0.3044 


0.7102 


0.008 


6.65 


410" 


2.2603 


0.8811 


0.2927 


0.0022 


0.0790 


0.9855 


0.3340 


0.7794 


0.008 


6.68 


410" 


2.4881 


0.4758 


0.1581 


0.0011 


0.0813 


1.1627 


0.3940 


0.9343 


0.008 


6.70 


410" 


2.4949 


0.2760 


0.0917 


0.0000 


0.0776 


1.2372 


0.4192 


1.0078 


0.008 


6.72 


410" 


2.3095 


0.1077 


0.0358 


0.0000 


0.0718 


1.2473 


0.4227 


1.0720 


[SI I] 


[SI 1 1] 


[SI 1 1] 


[SIII] 


[Nell] 


He I 


He I 


Hell 


Ha 


log hHp 


logw 


6731 A 


6312A 


9069A 


9532A 


3869A 


4471 A 


5836A 


4686A 


6563A 


ergs" ' 





0.0995 


0.0140 


0.2966 


0.7357 


0.2600 


0.0336 


0.0843 


0.0000 


2.9635 


38.727 


-0.56 


11.00 


0.1050 


0.0149 


0.3013 


0.7472 


0.2848 


0.0335 


0.0841 


0.0000 


2.9606 


38.734 


-1.03 


19.00 


0.1081 


0.0154 


0.3051 


0.7567 


0.2967 


0.0335 


0.0840 


0.0000 


2.9590 


38.741 


-1.21 


23.64 


0.1117 


0.0161 


0.3128 


0.7756 


0.3115 


0.0335 


0.0839 


0.0000 


2.9566 


38.753 


-1.37 


28.69 


0.1121 


0.0157 


0.3170 


0.7862 
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where e is the total injected mechanical energy (SN and stel- 
lar winds) per unit time in units of 10^' ergs s"' , ;i is the interstellar 
medium density in units of cm"', and t is the age of the shell in units 
of 10"* yr. Since we use the value of the energy at each time-step, 
this radius represents the instantaneous size of the region obtained 
by adding the energy produced by stellar winds and SN from one 
time step to the next. Therefore e is this energy divided by the time 
step. We have used the radius given for each time t and 1 Mq in Ta- 
ble 4 from Paper I, scaled to the stellar cluster mass of each model 
extrapolating this bubble geometry to a shell structure formed by 
the combined effects of the mechanical energy deposition from the 
massive stars winds and SN explosions belonging to the ionizing 
cluster. The ionized gas is assumed to be located in a thin spherical 
shell at distance Rj from the ionizing source. This approach has the 
advantage of eliminating the ionization parameter as a free parame- 
ter in the models, being now computed from the physical properties 
of the evolving young cluster. This radius is still large enough when 
compared to the shell thickness as to validate the approximation of 
a plane-parallel geometry. 

The hydrogen density has been considered constant through- 
out the nebula and an ionization-bounded geometry has been as- 
sumed, with the hydrogen density equal to the electron density 
for complete ionization (Case B of H recombination). We have 
computed models with two different values of n^'- 10 cm"' and 
100 cm"', in order to check the density effect on the emitted spec- 
trum. A density o f 10 cm"' is representative of small-medium iso- 
lated Hii regions jCastellanos et al. 2002a; Pe rez-Montero & Dfazl 
l2005h w hile 100 cm ' is m ore appropriate for modelling Hit 

S ;alaxies I Hagele et al , 20081) and large circum nuclear Hn regions 
Garcia- Vargas et al.L 1 19971 : iDfaz et al.L |2007|) . frecuently found 
around the nuclei of starbursts and AGNs. Although the constant 
density hypothesis is probably not realistic, it can be considered 
representative when the integrated spectrum of the nebula is ana- 
lyzed. 

Using the number of ionizing photons, Q(H), the nebula size, 
Rj, and the hydrogen density, we can obtain the ionization pa- 
rameter, u, given by: 

am 

" = n (2) 

where c is the speed of light. Under these assumptions, the parame- 
ter u recovers its physical meaning since it is directly derived from 
the physical cluster's properties. Hence, for a given IMF and n^/, 
each model is computed as a function of the cluster mass, age and 
metallicity. The values of u derived for our set of models are within 
the range of observed values for Hii regions. 

Since the current set of models have been computed for 
medium to large extragalactic Hii regions, whose measured red- 
dening is usually very low and often consistent with the galactic 
reddening, we have not included dust effects, as explained before, 
except for the abundance depletion. Only when these models are 
used for obscured compact Hii regions or for large star forming 
regions like starbursts galaxies, dust effects should be taken into 
account. 



3 RESULTS 

3.1 Emission lines. 

Once the photoionization code is applied, we obtain the emission 
line spectra of the associated Hn region produced by the ionizing 
cluster for every set of models. We have followed the evolution of 



the emitted spectrum during 5.2 Myr after the cluster formation. 
After this time, no emission line spectrum is produced due to the 
lack of ionizing photons. In the present study we have considered 
only the most intense emission lines in the optical spectrum which 
dominate the cooling for moderate metallicities. In higher (over so- 
lar) metallicity regions the cooling is shifted to the infrared. 

The models results, assuming a SAL2 IMF for all masses 
and metallicities, are summarized in Table [2] a) and b), for val- 
ues of o 10 cm"' and 100 cm"' respectively. Listed in columns 
I to 3 are: Z (metallicity), logarithm of the age in years, clus- 
ter mass (in Mq); in columns 4 to 20 we give the intensities 
relative to H^ of the following lines: [OII]/}3727, [OIII]/15007, 
[OIII]/14959, [OIII]/14363, [OI].i6300, [NII]/16584, [NII]i6548, 
[SII].i67I6, [SII].i6731, [SIII]A6312, [SIII]i9069, [SIII].i9532, 
[NeIII]i3869, HeU4471, HeU5836, Hem4686, H„; finally, in 
columns 21 to 23 we list the logarithm of the intensity of H^ in 
erg- s"' {logLHfS), the logarithm of the ionization parameter (log ;() 
and the ionized region radius in pc (i?,). 

Part of Table|2ji for Z=0.008 and a cluster mass of 4x lO^* M^is 
shown as an example. The complete tables are available in elec- 
tronic format. 



3.2 Comparison with previous models 

In order to evaluate the effect of the changes incorporated into the 
new models we hav e compared our present resu lts with the pre- 
vious ones given by iGarcia- Vargas et al] l ll995al) by computing a 
set of models using the same input abundanc es and the same as- 
sumptio ns for the ionization parameter as in iGarcia- Vargas et al.l 
( Il995hl) . In those models the ionization parameter was calculated 
with a fixed radius of log Rs(cm) = 20.84. The results are shown in 
Fig.[T] where the emission line ratios as a function of the ionization 
parameter log m, for Z = 0.02 at ages of 2 Myr, 3 Myr, 4 Myr and 
5 Myr can be seen. We have plotted the results for the following 
emission fines: [0II]/13727, [OIII].i5007, [OI].i6300, [NII].l6584, 
[SII].i6731 and [SIII].i9069 . In these figures, th e dotted lines cor- 
respond to the models from I Garcia- Vargas et alj ( Il995bl) while the 
solid lines correspond to the models presented here. 

We can see tha t , in gen eral, the emission lines resulting from 
iGarcia- Vargas et alj ( Il995bl) models are more intense than in the 
models computed in the present work for the same value of the ion- 
ization parameter. This is due to the use in the new set of models of 
NLTE blanketed atmosphere models for massive stars, which pro- 
duce fewer hard ionizing photons. This new set of models explains 
in a natural way the emission line ratios found in low excitation 
high metallicity Hn regions, as we will show in Section 4 and there- 
fore do not require ad hoc explanations to keep the number of hard 
ionizing photons to the observed values. 



3.3 Emission-line evolution 

In this work we follow the evolution of a given cluster during 5.2 
Myr from their formation. After this time, in most cases and de- 
pending on metallicity, the gaseous emission lines are too weak to 
be measured due to the paucity of ionizing photons. 

Figs [2^) to |2^) show the time evolution of the H„ lu- 
minosity and several emission line ratios, as labelled, se- 
lected to represent the most relevant ionization stages of the 
most common elements: [OIII].Li5007,4959/H/j, [OII]/13727/H/j, 
[OII].^3727/[OIII].^.^5007,4959, [OI]/16300/H„ and [NII]^6584/H„ 
for models of different metallicity, computed with SAL2 IMF. 
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Figure 1. Evolution of emi ssion lines at Z = 0.02 and SALl IMF for nn 
iGarcia- Vargas et al.l fT995bh while solid lines correspond to new models. 
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Figure 2. a). Evolution of the H„ luminosity and five emission line ratios for Z = 0.0001 (SAL2). The ratios are [OIII]/ii5007,4959/H/3, [OII]/13727/H/j, 
[OII]/[OIII], [OI]/16300/H„ and [NII]/16584/H„. The more appropriate hydrogen recombination line has been chosen for normalization in order to minimize 
reddening effects. Solid lines correspond to models with n// = 10 cm"^ while dotted Unes correspond to models with n^ = 100 cm"'. Black, red, green, blue, 
yellow, brown and grey lines correspond to different cluster masses: 0.12, 0.20, 0.40, 0.60, 1.00, 1.50 and 2.00 X lO' MO, respectively. 



In each figure, different colours correspond to a difi'erent cluster 
masses, from 0.12 x 10^ to 2 x 10^ Mq. The two different density 
cases, n//=10 cm"'' and nH = 100 cm"', are represented by solid and 
dashed lines, respectively. 

For the lowest abundance models (Z< 0.0001) emission line 



ratios involving oxygen lines are very weak (note the different scale 
in each metallicity panel) and its evolution is smooth. For interme- 
diate metallicity cases, the emission lines are intense, increasing 
with cluster mass. 

In the first Myr of the evolution, the [Oni]]ii5007,4959 lines 
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Figure 2. b). Evolution of emission lines for Z = 0.0004 (SAL2). Different colors and line code have the same meaning as in Fig.[2ji. 
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Figure 2. c). Evolution of emission Unes for Z = 0.004 (SAL2). Different colors and line code have the same meaning as in Fig.|2^. 



are intense and then they decrease to rise again at 3-4 Myr, due to 
the presence of WR stars, that makes the equivalent effective tem- 
perature of the cluster become higher and causes an increase of the 
excitation during a short time. The appearance of the WR features 
and their duration depends on the mass and metalhcity of the cluster 
as discussed in Paper I. The [OII]/^3727 line shows moderate vari- 
ations with metallicity. The [OII]/H^ ratio increases after 3 Myr of 
the star formation, being almost constant for low metallicity mod- 
els, before this age. As a result, the [OII]/[OIII] ratio is very low for 



all metallicities at ages younger than 4 Myr. The oxygen emission 
lines are most intense for models with Z=0.004. We will come back 
to the oxygen line evolution in the next section, when studying the 
diagnostic diagrams. 

The low ionization line [OI]/i6300 is an indicator of the 
prese nce of strong winds and supemovae jStasinska & Leithereii 
Il996l) . and therefore, a burst age indicator. On the sole basis on the 
photo-ionization models, the intensity of this line increases with 
cluster age and [0/H] abundance, therefore it shows higher val- 
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Figure 2. e). Evolution of emission lines for Z = 0.02 (SAL2). Dilferent colors and line code have the same meaning as in Fig.|2t. 



ues in old, evolved and high metallicity Hn regions. The addition 
of shock contributions (not included in our models), would render 
this value even higher. 

Finally, the [NII]/16584/Hq' line intensity ratio, which is in 
principle a good metallicity indicator, increase with increasing 
metallicity, even for regions with metallicities higher than Z=0.004. 

We have also produced photoionization models for ages older 
than 5.2 Myr (until 10 Myr) for the lowest metallicity cases, 



Z=0.0001, 0.0004 (see Fig.[3j. The lower the metallicity, the hot- 
ter the ZAMS and the longer the time spent by stars in the ZAMS. 
Hence, for the lowest metallicity clusters (Z=O.OOOI, 0.0004) we 
have explored the evolution further, to 10 Myr. However, the ioniza- 
tion parameters obtained are low, producing very weak (and some- 
times undetectable) emission lines from high-ionization stages, like 
[OIII], specially for the lowest masses. Of course, these regions do 
exist and produce a low-ionization spectrum, with observable [OI] 
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Figure 3. The ionization parameter as a function of age (beyond 5.2 Myr) 
for a 1 X 10^ Mo cluster, low metallicities (Z = 0.0001 and Z = 0.0004), 
and uh = 10 cm"' 




1 (Myr) t (Myr) 

Figure 4. Extendended evolution (uo to 10 Myr) of the H(, luminos- 
ity and five emission line ratios: [OIII]]/li5007,4959/H;3, [OII]/l3727/H/j, 
[On]/[OIII], [OI]^6300/H„ and [NII]^6584/H„, for a 1 X 10^ Mq cluster, 
with Z=0.0001 (SAL2) and n^ = 10 cm^^ 




Figure 5. Extendended evolution (uo to 10 Myr) of the Hn. luminos- 
ity and five emission line ratios: [Om]]/li5007,4959/%, [OII]/!3727/H/j, 
[OII]/[OIII], [OI]^6300/H„ and [NII]^6584/H„, for a 1 X 10^ Mq cluster, 
with Z=0.0004 (SAL2) and n« = 10 cm"'. 
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Figure 6. The [SII]/i/i6717,31/H„ line ratio vs age (in Myr) for low metal- 
licities (Z=0.0001 and Z=0.0004) and all masses. Solid lines join models 
with Z=0.0001 while dotted lines join models whith Z=0.0004. Dift'erent 
colors correspond to different values of the cluster mass. Models with val- 
ues of uh = 10 cm"' have been plotted in the upper panel while models 
with nn = 100 cm"' have been plotted in the lower one. 

and [Oil] and [Nil] lines but negligigle [OIII] line intensities which 
place them out of the plotted area of the diagnostic diagrams used 
in this work. Figures|4]and[5]illustrate this fact showing the results 
of the photoionization models for the 10 Myr evolution of a ioniz- 
ing cluster of 1 X IO^Mq and a density of = 10 cm"' for two 
metallicities: Z=0.0001 and Z=0.0004. 

Figs. |6] and |7] show the resulting sulphur emission line ra- 
tios: [SII].1.16717,31/H„ and [SII]/l/16717,31/[SIII]/l/19069,9532, as 
a function of the burst age. Sulph ur line ratios were proposed as 
ionization parameter indicator by lOfaz et al] l ll99lh and calibra- 
tions based on photoionization models at different metallicities 
were presented in ??. We have carried out a recalibration of this 
parameter using the models presented here that is detailed in Ap- 
pendix A. Since the ionization parameter is directly related to the 
stellar cluster age, as suggested by Garcia- Vargas et al. ( 19953), we 
propose to use these sulphur lines as age indicators for Hii regions. 
Fig. [6] shows the [SII]/i/i6717,31/H„ ratio for the lowest metallici- 
ties (Z=0.0001 and Z=0.0004) and all masses. Models with differ- 
ent values of Yin have been plotted in different panels. For a given 
metallicity, the plot shows that the line ratios involving the sulphur 
lines are low and rather constant during the first 3 Myr of evolution. 
From 3 to 5 Myr, they increase smoothly with age for all cluster 
masses. The line ratios also increase with cluster mass at a given 
age but cluster mass can be constrained via the derived number of 
ionizing photons from B aimer line luminosities. 

3.4 Emission line intensities and the ionization parameter. 

Figs. [8] a) to e) show the intensity, with respect to of the 
following emission lines: [OII]-i3727, [OIII].15007 , [OI].16300, 
[NII]i6584, [SII]^6731, [SIII]^9069 vs logu, obtained from the 
computed models, for metallicities Z = 0.0001, 0.0004, 0.004, 
0.008 and 0.02. In each figure four panels are shown for ages of 2 
Myr (top left), 3 Myr (top right), 4 Myr (bottom left) and 5 Myr 
(bottom right). Solid lines correspond to = 10 cm"' models 
while dotted lines correspond to Uu = 100 cm"'. In total 7 points 
per line are shown corresponding to 7 values of the cluster mass: 
0.12, 0.20, 0.40, 0.60, 1.00, 1.50 and 2.00 xlO^ Mq. We can see 
in the figures that models with njj = 10 cm"' and 100 cm"' follow 
the same trend, and cover different ranges of the ionization param- 
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Figure 8. a). Evolution of emission line intensities, as labelled, noiTnalized to the Hp intemsity for Z = 0.0001 (SAL2) models as a function of the logarithm 
of the ionization parameter. Solid lines correspond to models with n// = 10 cm"^ and dotted Unes correspond to models with n// = 100 cm"'. 
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Figure 8. b). Same as Fig.[8]but for Z=0.0004. 



eter value, with this being higher, as expected, for the lower density 
case. 

For a given metallicity, the changes in the emission line spec- 
trum and in the ionization parameter are due to the changes in the 
cluster mass, which determines the number of ionizing photons, 
and in the cluster age, which influences not only the total number 
of ionizing photons but also the overall ionization spectrum hard- 
ness and the ionized region size. For example, the [OIII]/i5007 line 



is intense during the first few Myr of the cluster evolution, for all 
metallicities, and then it decreases, to rise again at 4 Myr due to the 
presence of WR stars that produces a harder ionizing continuum. 

From our models we predict an intrinsic sizing elfect in the Hii 
region evolutionary sequence, which implies a natural decreasing 
of the ionization parameter. For a very young region whose ioniza- 
tion is dominated by the 0-B stars, the size is still small but the 
number of ionizing photons is high, resulting in a high ionization 
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parameter. As the cluster evolves, the region size increases due to 
the stellar winds, but the number of ionizing photons decreases only 
slightly producing a decreasing ionization parameter, and there- 
fore lower intensity lines. For the highest metallicities, this effect is 
more remarkable. High metallicity regions have less ionizing pho- 
tons and more intense winds therefore implying a lower ionization 
parameter. In addition, the effective cooling of the gas is transferred 
from the optical emission lines to the near and mid infrared ones as 



metallicity increases. These two facts explain the fact that in high 
metallicity regions, the optical emission lines are extremely weak, 
reaching in most cases only about one per cent of the Yip line inten- 
sity. These lines are difficult to detect and this fact has biased the 
observational samples during a long time towards regions of a re- 
stricted metallicity range which has important implications for the 
interpretation of photometric samples of Hii regions. 
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Figure 8. e). Same as Fig.[8]but for Z=0.02. 




Figure 7. The [SII]/l.i6717,31/[SIII]/U9069,9532 line ratio vs age (in Myr) 
for iiigii metallicities (Z=0.004, 0.008. 0.02) and all masses. Solid lines join 
models with Z=0.004, dotted lines those with Z=0.008 and dashed lines 
those with Z=0.02. Different colors correspond to different values of the 
cluster mass. Models with values of n/y = 10 cm"^^ have been plotted in 
the upper panel while models with n^ - "^'^ nn,-^ 
lower one. 



100 cm have been plotted in the 



4 DISCUSSION 

4.1 Observational Data 

We have used first a compilation of data on Hii galaxies, high and 
low metallicity Hii regions, and circumnuclear star forming regions 
regions (CNSFR), with emission fine intensities measured both for 
the auroral and nebular [SIII] lines at /i 6312 and AA 9069,9532 re - 
spectively. The data for CNSFR can be found in iDiaz et al.l jioOTh 
where the appropriate references for the rest of the objects in the 
compilation can also be found. For all the regions, the electro n den- 
sity has been derived using the [SII] ratio ( IOsterbroc3 . 1 1 9891) . 

For this sample, metallicities have been derived following 



standard techniques in the cases in which electron temperatures 
could be derived. Othe rwise, an empiric al calibration based on the 
S23/O23 parameter (see lDiaz"et has been used. 

Data from ICastellanos et all ( l2002al) correspond to Hii re- 
gions in the spiral galaxies: NGC628, NGC925, NGC1232 and 
NGC1637. The Hii regions have been splitted by density, which 
was derived from the [SII] ratio. For seven regions, ion weighted 
temperatures from optical forbidden auroral to nebular line ratios 
were obtained, and for six of them oxygen abundances were de- 
rived using empirical calibration methods. For the rest of the re- 

i ;ions, metallicities have be en estimated from the S23 calibration 
Perez-Montero et al. and tailored photoionization models. 

iGarcia- Vargas et al.l [l99% give data of 4 Hii giant circum- 
nuclear regions of NGC 7714. As usual, densities were derived 
from the [SII] lines. Oxygen abundances were obtained by the 
electron temperature method thanks to the direct detection of the 
[OIIIj/14363 Une. 

IZaritski. Kennicutt & Huchr3 ( Il994h made an analysis of 159 
Hii regions in 14 spiral galaxies, from which we have chosen those 
that have well measured [SII]]^^6717,31 and [SIII]]/1/19069,9532 
emission lines to obtain th e density and the oxygen abundance 
through the S23 parameter jPerez-Montero et"aL i l2006h . This re- 
duces the number of HII regions to 36. 

We have also taken data from Ivan Zee & Havnej l [2006h , cor- 
responding to 67 Hii regions in 21 dirr. They provide the emis- 
sion line intensities of the [0II]/13727, [OIII]/l/l5007,4959 and 
[SII]/l/l6717,31 lines. In addition they also provide the [SII] line 
ratio, which we have used to derive the electron density. The [OIII] 
lines have been used to obtain the electronic temperature and con- 
sequently the oxygen abundance. 

Data from llzotov & Thuanl ( |2004|) consist of Hii regions in 76 
blue compact dwarfs (BCD) galaxies whose abundances have been 
derived using electron temperatured from the ratio of the nebular 
to auroral lines [OIII] lines. All the regions of this sample have 
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Table 4. Abundance ranges for the comparison of models and observational 
data. 



Label 


Observed Metallicity Range 


Z in plotted models 


1 


12+log(0/H) < 7.7 


0.0004 


2 


7.7 < 12+log(0/H) < 8.4 


0.0040 


3 


8.4 < 12+log(0/H) < 8.7 


0.0080 


4 


8.7 < 12+log(0/H) 


0.0200 



oxygen abundances lower than 12 + log(0/H) = 8.5. Data have 
been separat ed by density by means of the [SII] line ratio. 

Finallv. lYin et al.l l l 2007h provide data for 53 1 galaxies and Hii 
regions from the SDSS-DR4 sample, from which the [OIII]/i4363 
line has been measured and the oxygen abundances have been de- 
rived by standard methods. These data do not include the [SII] lines 
and hence no electron density could be derived. The different ob- 
servational data sources and the method used to derive the oxygen 
abundances are summarized in Table [3] 



4.2 Diagnostic Diagrams: optical emission line ratios 

We have plotted together the model results and the observational 
data in different diagnostic diagrams in figures from|9]to[T2l These 
diagrams can be used to study the relationship among different 
emission line ratio and to extract information about the physical 
parameters of the ionizing clusters. 

Each figure is divided in 4 panels, one per metallicity: Z = 
0.0001 and 0.0004 (top left panel), 0.004 (top right panel), 0.008 
(bottom left panel) and 0.02 (bottom right panel). The highest 
metallicity value of the grid, Z = 0.05, has been excluded due to 
the absence of spectroscopic data with those high abundances. The 
lowest metallicity model tracks are plotted in the same panel as 
those corresponding to Z=0.0004 to compare the position in the 
diagram of the emission lines of such low metallicities. Solid and 
dotted lines represent models with n//=10 cm"^ and n/;=IOO cm"' 
respectively. In each panel the model metallicity is the central value 
of the different ranges, except for ranges 1 and 4 that include data 
with metallicities lower than Z = 0.0004 and higher than Z = 0.02 
respectively. Different line tracks correspond to the evolution of 
the models with different cluster masses, from the lowest one on 
the right to the highest on the left. The observational data, taken 
from different sources as explained in the previous section, whose 
derived metallicities are in therange of the models we want to com- 
pare with, are shown in each panel. They have been divided by den- 
sity and subdivided by metallicity according to the ranges given in 
Table m 

The Hit regions plotted here describe a well defined se- 
quence in diagnostic diagrams which involve oxygen line ra- 
tios, like log([OIII]/H^) vs log([OII]AI/3) and log([OIII]/H/i) vs 
log([OII]/[OIII]) ( Fig. |9] and Fig. [TUll. Most of the observed re- 
gions are located in the panels corresponding to Z=0.004 and 
Z=0.008 since few data exist for very low or oversolar metallici- 
ties. In fact, this is a bias effect since most popular Hit regions have 
been detected and catalogued through their high intensities of the 
[OIII]/l/l5007,4959 lines, which happen in medium-low metallic- 
ity. Furthermore, it is at these metallicities, Z=0.004 and Z=0.008, 
that a direct measurement of the [OIII]/i4363 line, used to derive 
the oxygen abundance, is possible. 

Models of Z = 0.0001, 0.0004 and Z = 0.02 do not have the 
aim of reproducing exactly the observation sequence but should be 
understood as metallicity limits for the observed regions. Observa- 



tional data in these panels correspond to regions with metallicities 
lower than 12+log(0/H) = 7.7, in the first panel (upper-left), or 
higher than 12-l-log(0/H) = 8.7 in the last panel (bottom-right), as 
summarized in Table |4] 

There are no observational data with Z < 0.0004. These re- 
gions have very low values of the [OIII]/i4363 line, undetectable 
with moderate size telescopes and standard S/N values, and there- 
fore their abundances have not been derived. 

Our models predict that very low metallicity regions have 
small sizes in Ha images where most of these Fin regions are cata- 
logued and therefore observational samples can be seriously biased 
towards medium to large size Hn regions with metallicities Z > 
0.001 (see Section 4.4 and fig. [14] for a detailed discussion). For 
high metallicities, the [OIII]/i4363 line is not strong enough to be 
detected at a level of T^ determination and abundances have been 
derived through alternative methods as already described. 

Our models cover the observational data range corresponding 
to different cluster masses. Different masses imply different num- 
bers of ionizing photons and hence different values of the ionization 
parameter. These values decrease with time as the clusters evolve. 
The models involving the most massive stellar clusters have the 
highest ionization parameters, due to their higher number of ioniz- 
ing photons at a given age. Less massive clusters follow a similar 
time evolution to the most massive ones, but starting at a lower ini- 
tial value. 

As explained before, regions with intense 
[OIII]/l/15007,4959/H,j line ratios correspond to the highest 
ionization stages. For a given cluster mass, as the cluster evolves, 
the high excitation emission lines tend to decrease while [OII]/H^ 
increases. This is due to the change of the ionization parameter 
with cluster evolution. For a very young region, whose ionization 
is dominated by 0-B stars which produce a great amount of 
ionizing photons , the size is still small, and the resulting ionization 
parameter is high. As the cluster evolves, the region size increases 
due to the stellar winds while the number of ionizing photons 
remains almost constant (or even decreases) with respect to early 
phases. This effect makes the ionization parameter, hence the 
emission line intensities and in particular the oxygen lines, to 
decrease. 

This evolution influences also the data sample since early ages 
result in more compact regions (therefore high surface brightness) 
with high ionization parameters. Thus, the emission line intensi- 
ties of the [OIII]/l/l5007,4959 lines are higher, increasing their de- 
tectability in a given sample obtained with the same instrumental 
configuration. In fact, most of the derived ages for those regions 
would be younger than 5 Myr, when they emit enough ionizing pho- 
tons to produce a detectable H„ luminosity and from those, most of 
them would be biased towrads younger ages (less than 3.5 Myr). 

In all Figs.|9]to [T3] there are not observations found in the tail 
of our models with high [OII]i3727 and low [OIII]/l^5007,4959 
intensities. These values imply very low ionization parameters, 
that indicates more diffuse regions with low excitation. Since in 
most samples only the brightest Hii regions have been selected, 
due to their conspicuous Ho' emission or their high values of the 
[OIII]/i/i5007,4959 line intensities the lowest excitation regions 
tend to be excluded. Some of these regions could also be older low 
metallicity regions (in the age range from 5 to 10 Myr, see Figs. [4] 
and |5j. 

The influence of the electron density, although of second or- 
der, should also be considered. The models with n//=I00 cm"' 
show lower values of the ionization parameter. However, the ra- 
diation has more particles per cm"' to ionize, and the excitation 
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Figure 9. Diagnostic diagram sliowing log([OIII]/i/i5007,4959, /H^) vs Iog([OII]/l3727/H^) for models computed witli SAL2 IMF. Each panel show s models 
for a d ifferent metallicity and the corre sponding observational data in tha t metallicity range. Observational data come from: |Zaritski. Kennicutt & Huchral 
h994 (circles), llzotov &Thuai]j2004l) (sq uares). ICastellanos et ai)2002bl) (diamonds). lvan Zee & Haviiesl i2006l) (triangles upi lbarcfa- Vargas et aljj 19971) 
(red triangles down), and IPiaz et^ (asteriscsof different colours: black: CNSFR; green: high metallicity HII regions; brown: low metallicity Hn 

regions; magenta: HII galaxies). Filled symbols represent objects with nn = 100cm"', while open symbols have nn = 10cm"', for data whose density has 
been obtained from the [SII] ratio. 
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Figure 10. Diagnostic diagram showing log([OIII]^5007.4959/H,j) vs log([OII]/i3727/[OIII]^,^5007,4959) for models computed with SAL2 IMF. Each panel 
shows a different metallicity model and the corresponding data in that metallicity range. Symbols have the same meaning as in Fig. [9] 
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Figure 11. Diagnostic diagram showing log([OIII]/{i5007,4959/H/j) vs log([NII]/l6584/[OII]/i3727) for models computed with SAL2 IMF. Each panel shows 
a different metallicity model and the corre sponding d ata in that metallicity range. Symbols have the same meaning as in Fig.|9] with the addition of the blue 
dots, which correspond to Hu regions from lYin et al] i 2007.) . 




Figure 12. Diagnostic diagram showing log([OIII]/l/i5007,4959/H/j) vs log([NII]/16584/Hj,) for models computed with SAL2 IMF. Each shell shows a different 
metallicity model and the corresponding data in tha t metall icity range. Symbols have the same meaning as the same than in Fig. |9] with the addition of blue 
dots which correspond to Hii regions from lYin et al ] j2007h . 
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Table 3. Observational data used for this work 



Reference 



Objects 



[O/H] method 



^fazetaL Q007, and references therein) 
^astellanoSj^^iaz^Terievich (2002a) 
^arcia;^^as_e£ah Q99T) 
Zarit ski. Kennicutt & Huchra (1994) 
^an Zee & Haynes (2 006) 
Izotov & Thuan (2004) 
Yin et al. ( 2007) 



High and low metallicity regions, CNSFR and Hii galaxies 
high metallicity Hii regions 
Hi: circumnuclear regions 
Hii regions in Spirals 
Hi: regions in dlrr 
Hii regions in BCD 
Hii regions 



Direct method and S23/O23 
Direct method and S23 
Direct T,, method 

S23 

Direct method 
Direct method 
Direct T,, method 




Figure 13. Diagnostic diagram showing log([SlI].i6717,6731/H„) vs log([OlIl]^/15007,4959/H;?) for models computed with SAL2 IMF. Each shell shows a 
different metallicity model and the corresponding data in that metallicity range . Symbols have the same meaning as in Fig.|9] 



degree is higher, producing model sequences shifted upwards from 
the lower density ones. More significant differences between mod- 
els of different densities are found in the highest metallicity dia- 
grams (Z=0.02). 



4.3 Sulphur and nitrogen emission line ratios. 

Diagnostic diagrams involving log([NII]/16584/H„) and 
log([SII]/i/i6717,31/H„) ratios have the advantage of using lines 
very close in wavelength, and therefore virtually reddening-free. 

In large star forming regions, like starburts galaxies, the rela- 
tionship between the excitation parameter [OIII]/1/15007,4959/HjB 
and the metallicity indicator [NII]/16584/H„ gives us informa- 
tion about the youngest populations plu s the star formation rate 
jKennicutt. Tamblvn. & Congdonl . [l994h . At the beginning of the 
evolution H„ emission is intense and comes from the ionizing pho- 
tons produced by the most massive stars, most of them still in the 
main sequence. This number of ionizing photons increases with the 
cluster mass and decreases with metallicity. As the galaxy evolves, 
the Hj, emission line decreases, while the [NII]/i6584 line increases 
because of the gas ejection by stars. 

The [NII]/t6584/[OII]/13727 ratio was introduced as a good 



metallicity indicator by iDopita et alj bOOOh . Nitrogen is a sec- 
ondary element, or, at least, it has a strong secondary component, 
which makes the ratio [NII]/[OII] to increase with metallicity. Be- 
sides, the value of the mean temperature of Hn regions decreases 
as metallicity increases. When metallicity rises, the electron tem- 
perature becomes too low to excite the transitions, hence [Oil] 
decreases while [Nil] does not, thus increasing the [NII]/[OII] ra- 
tio. A good correlation has been found between [NII]/[OII] and the 
N^/O ^ ionic abundance ratio , whic h is as sumed to trace the N/0 
ratio jPerez-Montero & Diazll2005l) . Also. lDi'az et al.l j2007h have 
shown that a very tight correlation between [NII]/[OII] and N/0 
exists for high metallicity Hn regions and circumnuclear star form- 
ing regions, which follow a sequence of increasing N/O ratio with 
increasing emission line ratio. 

In our diagnostic diagrams of Fig. [T2] we can see how this 
ratio changes with metallicity. For low metallicity models, both 
the [OIII]/H^ ratio and [NII]/[OII] ratio are low. As metallicity in- 
creases both line ratios also increase. When metallicity reaches the 
limit of Z=0.004, [OIII] starts to decrease while [NII]/[OII] ratio is 
still growing. At Z=0.02 the models show low values of [OIII] and 
high values of [Nil]. The only available observational data in this 
part of the plot are high metallicity Hn regions and CNSFR, show- 
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ing high values of [NII]/[OII], which in dicate high N/0 ab undance 
ratios, as can be seen in the sequence of lOfaz et al.l ( l2007h . 

Regarding the nitrogen abundance, it is well known that N is 
mainly primary for low metallicities and it is always possible to 
scale the N abundance that we include as input in cloudy by using 
a given function, such as the one obtained by iMolla et al] l l2006h 
for the relationship of N/0 vs 0/H. However, we remark that in the 
models presented in this paper, we use SSP models where chem- 
ical evolution has not been taken into account, either for N or for 
the other elements as Fe or Ca. The nitrogen abundance has been 
scaled to the solar value as it has been done for the rest of the ele- 
ments. Each element has its own evolution and in order to take this 
into account it is necessary to make precise calculations follow- 
ing a c onsistent method. Thi s method can be seen in our p revious 
works jMartm-Manion et al.l |2008| ; iMartin-Manioiil . |2009|) where 
we preent photoionization models that ise as imput abundances the 
results obtained from chemical evolution models. 

Finally, Fig. [13] shows the [SI1]/1/16717,31/H„ ratio vs 
[OIII]/(/(5007,4959/H/3 also almopst reddening-independent. We 
can see that high values of the [SII]/1/167 17,31 line imply a low 
ionization parameter, and therefore a low [OIII] emission line in- 
tensity. However, the [OIII]/1/15007,4959 lines are very metallic- 
ity dependent, and a low intensity of this emission line could be 
due to a high metallicity rather than to a low ionization parameter 
produced by an evolved population. Although both ratios are in- 
dicators of the ionization parameter, the values of ii derived from 
oxygen emission line ratios can be lower than ii the values derived 
from sulphur emission line ratios, which are not so metallicity de- 
pendent. 
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Figure 14. Evolution of the [OIII]/14363 /Hp ratio for each metallicity, clus- 
ter masses from 0.12 to Ix 10^ Mq and nH = 10. 
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4.4 Abundance determinations in Hn regions and model 
consistency 

Emission lines allow the determination of abundances by a direct 
method based on the electronic temperature, Tg, derived from suit- 
able emission lines from which those of [0111] are the most popular 
We have obtained oxygen abundances from the modelled emission 
lines in order to check which is the most effective method to recover 
the oxygen abundance used as input for the models. An accurate 
analysis requires the detection of the auroral line [01I1]/14363 with 
an intensity larger than 5 per cent of H/j. 

Fig. [T4] shows the evolution of the [OIlI]/14363 /H/j ratio for 
each metallicity, n//=10, and for cluster masses of 0.12, 0.20,0.4, 
0.6 and Ix 10^^ Mq. This ratio is high for Z=0.0004 and Z=0.004, 
close to 5 per cent of H^ in the first Myr of the evolution. However, 
for both, Z=0.0001 and Z=0.008, [OIII]i4363 /Hp has comparable 
values and below standard detection limits. Thus, the non detection 
of this line can indicate either a high metallicity or a low oxygen 
content. For Z=0.02 there is no [OIIIJ/14363 line produced for this 
density case. 

Fig. [15] shows the comparison between model oxygen abun- 
dances and the ones obtained by empirical methods. It can be seen 
that, from our photoionization models, we obtain the auroral line 
[01I1]/14363 along the whole cluster evolution for metallicities in 
the range 0.0001< Z < 0.008. The [OIII] line electron temperature 
method gives good results for all of them, with differences lower 
than 0.05 dex between input and derived abundances. 

At Z = 0.008, corresponding to an oxygen abundance of 
12+log(0/H) =8.5, the [OIII],i4363 line intensity range goes from 
1.62 per cent of the H^ intensity 1 Myr after the star formation, to 
0.82 per cent 4 Myr after the burst. However, despite this line be- 
ing detectable at this metallicity, the electronic temperature method 



5 5.5 6 6.5 

log t 

Figure 15. Comparison between abundances used as input for the models 
(solid lines) of different metallicities and abundances obtained using differ- 
ent empirical calibrators (dashed lines). 

is not able to reproduce the input metallicity of the model. Simple 
photoionization models are not able to reproduce the derived T^ in 
some situations. ^Stasinska (_1978» 2 002 . 2005j) showed that for high 
metallicity Hn region models, the oxygen abundance derived from 
predicted emission line ratios and the T^-method is lower than the 
one used as input in the models. These differences may be signif- 
icant when large temperature gradients or fluctuations are present. 
In these cases, the temperature obtained from the emission line, 
T;[OIII], is different from that derived from the [OIII] line ratio, 
T„ and from the average ionic temperature T(0^^), with T;[OIII] 
and T,- obtained from the models not being the representative tem- 
perature of the high ionization O ^^ zone. Fig.|16|shows a plot sim- 
ilar to Fig.l of IStasifiskal l l2005h comparing a sequence of model 
derived oxygen abundances as a function of the input values . In 
this diagram, we can see that both sets of values coincide for low 
metallicities; for metallicities around 8.6 however, deviations be- 
come notorious and for higher metallicities the deviations are much 
more significant. 

The abundances analysis in high metallicity Hit regions is 
hampered by the fact that oxygen optical lines act as main coolant 
for the nebula, hence a higher oxygen abundance leads to a more 
effective cooling. As the gas cools down, the electron temperature 
gets lower and the optical [OIII] forbidden lines get weaker. In gen- 
eral, this low excitation makes any temperature-sensitive lines in 
the optical range too weak to be measured. For this reason empiri- 
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Figure 16. Comparison between abundances used as input for the models 
and abundances obtained using T(,-based methods. Each color corresponds 
to a different age sequence. The values of 12+log(0/H) derived empirically 
are very close to the input values up to 8.6. 

cal calibrators i nvolving other lines besides the optical oxygen o nes 
should be used jPerez-Montero & Diazll2005l : lDiaz et alil2007l) . In 
our models, for metallicities Z = 0.008 and Z = 0.02, abundances 
have been determined with the ratio S23/O23, obtaining reasonably 
good results for Z=0.008 (with a difference of 0.2 dex approxi- 
mately), but not for Z=0.02, which shows important differences at 
the beginning of the evolution. For Z = 0.05 we cannot recover 
the oxygen abundance used as input by any empirical method since 
an effective calibration for such high metalli city range (which is 
appro ximately 3Zo in the reference system of iGrevesse & Sauvall 
Il998l) has not been yet devised. 

4.5 Ionized Region equilibrium time and ionizing cluster age 

The output of the photoionization models shows the picture of the 
ionized region once the equilibrium state is reached. This means 
that, when we observe a Hit region, what we really see is the effect 
of an ionizing cluster of a given age, regardless the actual age of 
the ionized nebula. The time step in our models corresponds to the 
stellar cluster age step given by the isochrone time resolution. The 
consequence is that the age resolution with which we can date the 
ionizing cluster from the information of the gas ionizing spectrum 
depends on the speed at which the nebula reaches the equilibrium, 
and this process is controlled by the metallicity. In the worst case, 
if the equilibrium time is longer than the cluster's age, it would 
be possible to see the effect of the ionizing cluster several Myrs 
after the ionization. Therefore, the ionizing stellar cluster would 
be older than inferred from the emission line spectrum, being this 
observed cluster the result of its evolution during equilibrium time 
of the nebula. [Wofford ( 2009;) shows that the time it takes for a Hit 
region to reach thermal equilibrium at a metallicity of Z = 0.001 
can be longer than 1 Myr within the first 5 Myr of evolution. The 
equilibrium time scale increases with decreasing hydrogen density. 
In addition, it increases at lower metallicities, due to the low content 
of metals for cooling the gas. 

We have tested this effect in our models by checking the equi- 
librium time in each cloudy model. Fig.[l7]plots the age of the ion- 
izing region given by cloudy once the equilibrium has been reached 
versus the age of the ionizing star cluster. Lines in orange, black. 
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Figure 17. Relationship between the age of the ionizing stellar cluster and 
the recombination age of the ionized region. Orange, black, red and blue 
lines join models for Z = 0.0001, 0.0004, 0.004 and 0.02 respectively. Solid 
lines con'espond to models for a cluster mass of to Ix lO' M© while dotted 
lines correspon to a cluster mass of 0.12x10^ Mo 



red and blue correspond to models with Z = 0.0001, 0.0004, 0.004 
and 0.02 respectively. Solid lines refer to models for a cluster mass 
of Ix 10'' Mo while dotted lines to models for a cluster mass of 
0. 12x 10^ Mo. All models have been computed with a density of 
n;y = 10 cm"\ 

A first conclusion to be extracted from this graph is that the 
age of the ionized region does neither depend on the cluster mass 
nor on the cluster age but it is controlled by the metallicity. This 
implies that neither the total number of ionizing photons nor the 
ionizing spectrum shape have a significant role in setting the equi- 
librium time of the ionized regions, which is controlled by the elec- 
tron temperature (given by the metallicity). This is clearly seen in 
Fig. [17] where dotted and solid lines overlap and the ionized region 
age is almost constant along the cluster evollution. According to the 
computed models, the average equilibrium age is reached around 
t,,= 13.6, 3.1, 0.2, 0.1 and 0.05 Myr for metallicities Z = 0.0001, 
0.0004, 0.004, 0.008 and 0.02, respectively. 

The second conclusion is that there are two metallicity 
regimes in which the ionizing cluster age has to be interpreted in a 
different way from the nebular emission line spectrum. In moder- 
ate to high metallicity Hn regions (Z > 0.004) the ionized nebula 
equilibrium time is shorter than the ionizing cluster age by a factor 
between 4 and 10. This means that for regions with metallicities 
above 0.004 we will derive the same age for the ionizing cluster 
whether we look directly to the embedded cluster (whenever pos- 
sible) or we observe the cluster effects on the surrounding nebula, 
that is the nebular emission line spectrum. In these cases, the res- 
olution time, and therefore the error with which we can determine 
this age, is around 0.05 Myr at Z=0.02, 0. 1 Myr at Z=0.008 and 0.2 
Myr at Z=0.004. In this metallicity regime, if we detect WR bumps 
in the optical spectrum or O, B and WR stars from UV images, we 
can assure, according to our results, that these stars are the ones 
responsible for the detected nebula ionization. On the contrary, in 
low metallicity Hn regions, there is a delay between the time at 
which we observe the ionized region and the age of the stellar clus- 
ter responsible for this ionization. At Z = 0.0004 this delay is of the 
or der of 3 Myr, w hat is consistent with the value of 1 Myr derived 
bv lWoffordl ( l2009l) for Z=0.001 . The delay can be as long as 13 Myr 
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Figure 18. Flow Diagram to unveil the ionizing cluster physical properties thi'ough the gas optical emission line spectrum 



at Z = 0.0001. This implies that we may observe low metallicity re- 
gions apparently ionized by a young cluster but without finding the 
expected young and hot stars when we observe directly inside the 
nebula. Moreover, we could find more evolved stars radiating in the 
IR as RSGs as the resulting products of the evolution of the ionizing 
stars. We suggest to take into account this result when interpreting 
the stellar content inside low metallicity Hii regions. 

4.6 Unveiling the physical properties of the ionizing clusters 

The purpose of this piece of work is to show that it is possible to 
unveil the properties of the young cluster responsible for the ion- 
ization on the basis of the emission line spectrum of the ionized gas 
alone. In what follows we suggest some recommendations to find 
the physical properties of the dominant young ionizing population 
(see Flow Diagram in Fig.llSt. 

- The first parameter than can be determined is the electron 
density, which is derived from the [SII]/l/l6716, 6731 lines. This 
value influences the ionization parameter and therefore it is an im- 
portant input for our models. 

- The metallicity can be derived or constrained once the elec- 
tron density has been determined. If the value of [OIII]/i4363 
/H^ is higher than 0.023, then abundances can be derived by 
standard methods using the electron temperature calculated from 
the [OIII]/l4363, [OIII]/l4959, [OIII].i5007 line ratio. We remark 
that the absence (or a very low value) of optical [OIII] lines is 
an indicator of very low or very high metallicity. PopStar in- 
cludes models with metallicity as low as Z=0.0001. If the value 
of [0III]/14363 fHp is lower than 0.023, we recommend to check 
the value of log([OII]/i3727 /H^). Values lower than 0.5 indicate 



very low metallicities while the opposite indicates that the region 
is in the high metallicity regime. In this latt er case, we can estimate 
the m etallicity through empirical methods IPerez-Montero & Dfazl 
l2005i) . These methods allow to estimate abundances up to solar val- 
ues, (o (Z = 0.02). Values higher than those are usually not found in 
optically selected samples since these regions are efficiently cooled 
by infrared lines and have a very faint optical emission line spec- 
trum. 

- The age should be constrained once the metallicity evalu- 
ation has been carried out. We propose to use the diagnostics 
with [SII]^^6716,31/H„ (Fig.|6), for low metalhcity regions, and 
[SII]-i-i6716,3I/[SIII]^^6069,9532 (Fig.|7), for high metalhcity re- 
gions. In both cases, we need to use the electron density to select 
among different sets of models. We have computed models for two 
density values: 10cm"' and 100 cm"', which we recommend to use 
when electron density is less than 50cm"' or higher than 50cm"' 
respectively, and up to a maximum value of 400 cm"' . On the ba- 
sis of these diagrams and the observed values of the line ratios, we 
can distinguish if the ionizing cluster is younger or older than 3 
Myr. The resolution at which we can estimate the age depends on 
the equilibrium time as discussed in section 4.5. In the case of low 
metallicity regions (below Z = 0.004) the equilibrium time is of 
the order, or even longer than the ionizing cluster age. This has to 
be taken into account when observations of the embedded stars are 
available, since in this low metallicity scenario these stars are the 
result of the evolution of the ionizing cluster and not the ionizing 
cluster itself. 

- The luminosity of H„ emission line is directly related to the 
number of ionizing photons, and therefore to the mass of the clus- 
ter. However, H„ luminosity depends not only on the cluster mass 
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but also on the age and metallicity, even during the first 5 Myr of 
evolution (see Fig.|2]l when the emission line spectrum in the opti- 
cal is detectable, according to our models. Moreover H„ is detected 
up to an age of 20 Myr, when the emission line spectrum is not yet 
detectable (we remind the reader that this paper is focused on the 
first few Myr when emission lines are prominent since these are the 
regions available in most Hii region catalogues). Due to the former 
reasons, it is difficult to derive the mass without previous estimates 
of metallicity and age. Therefore we propose to derive the mass 
of the cluster from the luminosity (Fig. [2} using the previously 
found restrictions on metallicity and age. 

- Finally, we can use the values of the complete set of emission 
lines to check the consistency of our determination. 

Since this process is slow to be applied to a large sample of 
objects, we are developing a fast software tool, based on the previ- 
ous guidelines and using minimization algorithms, to find the best 
model that fits the emission line values and determine the physical 
properties of the young cluster responsible for the ionization of a 
given Hit region, when spectroscopic data are available. 



5 SUMMARY AND CONCLUSIONS 

We present a new grid of evolutionary models as a tool to estimate 
the mass, age and metallicity of the dominant ionizing stellar popu- 
lation of a given Hn region from its optical emission line spectrum. 

We provide a complete set of photoionization models for Hii 
regions along the first 5.2 Myr of their evolution, using as ioniz- 
ing source a very young star cluster obtained from the theoretical 
evolutionary synthesis models PopStar (Mollaetal, 2009). |3 The 
basic grid is composed by SSPs computed with different IMFs and 
a wide range of metalficities (Z=0.0001, 0.0004, 0.004, 0.008, 0.02 
and 0.05). A Salpeter IMF with m,„„.=0.15Mo and m,,,, = IOOMq 
has been used in this work. The ionizing cluster is assumed to form 
in a single burst. Our grid considers five values of the total cluster 
mass (0.12, 0.20, 0.40, 0.60, 1.00, 1.50 and 2xlO'Mo), covering 
the range of observed Hii regions. As the cluster evolves, the stellar 
winds and supernovae blow the material, determining the Hii region 
size. This material is ionized by the emergent continua. We use the 
photoionization code cloudy to obtain the corresponding Hii region 
emission line spectrum once the equilibrium is reached. 

We present diagnostic diagrams where we compare our model 
results with a large sample of observed Hii regions, whose metal- 
licities have been consistently computed in a homogeneous man- 
ner according to the most appropriate method for each metallicity 
range. It is found that, in general, our models can satisfactorily re- 
produce the observations. 

Finally, we show that the emission line spectrum from an ion- 
ized region provides a powerful tool to unveil the physical proper- 
ties of the embedded young cluster responsible for the gas ioniza- 
tion and propose some guidelines to constrain its physical proper- 
ties. 

6 ACKNOWLEDGMENTS 

This work has been partially supported by DGICYT grant 
AYA2007-67965-C03-03 and partially funded by the 

^ More evolved cluster emission line spectra can be provided under re- 
quest, however, these spectra will have only low ionization emission lines 
and in particular do not present [OIII] lines. 



Spanish MEG under the Consolider-Ingenio 2010 Pro- 
gram grant CSD2006-00070: First Science with the GTC 
(http://www.iac.es/consolider-ingenio-gtc/). Also, partial support 
from the Comunidad de Madrid under grant S-0505/ESP/000237 
(ASTROCAM) is acknowledged. The authors thank an anony- 
mous referee for valuable comments and suggestions that have 
contributed significantly to the quality of this paper. 



REFERENCES 

Antia H. M., Basu S., 2005, ApJL, 620, L129 

Asplund M., Grevesse N., Sauval A. J., 2005, in ASP Gonf. Ser. 

336: Gosmic Abundances as Records of Stellar Evolution and 

Nucleosynthesis, 336, 25 
Bahcall J. N., SereneUi A. M., Basu S., 2005, ApJL, 621, L85 
Bressan A., Granato G. L., Silva, L., 1998, A&A, 332, 135 
Bressan A., Panuzzo P., Buson L., et al., 2006, ApJ, 639, L55 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Castellanos M., Di'az A. I., Terlevich, E., 2002, MNRAS, 329, 315 
Castellanos M., Di'az A. I., Tenorio-Tagle G., 2002, ApJ, 565, L79 
Castor J., McCray R., Weaver R., 1975, ApJ, 200, L107 
Clegg R. E. S., Middlemass D., 1987, MNRAS, 228, 759 
Clemens M. S., Bressan A., Panuzzo P., Rampazzo R., Silva L., 

Buson L., Granato G. L., 2009, MNRAS, 392, 982 
Chabrier G., 2003, ApJL, 586, L133 
Di'az A. I., Perez-Montero E., 2000, MNRAS, 312, 130 
Diaz A. I., Terlevich E., Pagel B. E. J., Vilchez J. M., Edmunds 

M. G., 1987, MNRAS, 226, 19 
Diaz A. I., Terlevich E., Vilchez J. M., Pagel B. E. J., Edmunds, 

M. G., 1991, MNRAS,253, 245 
Diaz A. I. 1994, in Violent Star Formation: from 30 Doradus 

to QSO's, Ed. G. tenorio-Tagle, (Cambridge University Press: 

Cambridge) p. 105 
Diaz A. I., Castellanos M., Terlevich E., Garcia- Vargas M.L., 

2000, MNRAS, 318, 462 
Di'az A. I., Terlevich E., Castellanos M., Hagele G. E, 2007, MN- 
RAS, 382, 251 
Dopita M. A., Evans I. N., 1986, ApJ, 307, 431 
Dopita M. A., Kewley L. J., Heisler C. A., Sutherland R. S., 2000, 

ApJ, 542, .224D, 
Dopita M. A., Groves B. A., Sutherland R. S., Binette L., Cecil 

G., 2002, ApJ, 572, 753 
Dopita M. A., Fischera J., Sutherland R. S., Kewley L. J., Lei- 

therer C, Tuff's R. J., Popescu C. C, van Breugel W., Groves 

B. A., 2006a, ApJ Supp. Series 167, 177 
Dopita M. A., Fischera, J., Sutherland R. S., Kewley, L., et al., 

2006b, ApJ, 647, 244 
Drake J. J., Testa P, 2005, Nature, 436, 525 
Ekstrom S., Meynet G., Chiappini C, Hirschi R., Maeder A., 

2008, A&A, 489, 685 
Evans I. N., Dopita M. A., 1985, ApJS, 58, 125 
Ferland G. J., Korista K. T, Vemer D. A., Ferguson J. W., King- 
don J. B., Verner E. M., 1998, PASR 1 10, 761 
Ferrini F, Penco U., Palla F, 1990, A&A, 231, 391 
Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 
Garcia Vargas M. L., Di'az A. I. 1994, ApJS, 91, 553 
Garcia- Vargas M.L., Bressan, A., Di'az, A. I., 1995a, A&AS, 1 12, 

13 

Garcia- Vargas,M.L., Bressan, A., Di'az, A. I., 1995b, A&AS, 1 12, 
35 



22 Martm-Manjon et al. 



Garcia- Vargas M. L., Gonzalez-Delgado R. M., Perez E, Alloin 

D. , Di'az A. I., Terlevich E., 1997, ApJ, 478, 112G 

Garcia- Vargas M. L., Molla M., Bressan, A., 1998, A&AS, 130, 
513 

Garcia- Vargas M. L., Molla M., Martm-Manjon M. L., 2009, MN- 

RAS, in preparation 
Garnett D. R., Dufour R. J., Peimbeit M., Torres-Peimbert S., 

Shields G. A., Skillman E. D., Terlevich E., Terlevich R. J., 1995, 

ApJ, 449, L77 
Grevesse N., Sauval A. J., 1998, SSRv, 85, 161G 
Groves B. A., Dopita M. A., Sutherland R. S., 2004, ApJS, 153, 

75 

Hagele G. F., Diaz A. I., Terlevich E., Terlevich R., Perez-Montero 

E. , Cardaci M. V., 2008, MNRAS, 383, 209 
ffiUier D. J., Miller D. L., 1998, ApJ, 496, 407 
Izotov Y. I., Thuan T. X., 2004, ApJ,602, 200 

Kennicutt Jr. R. C, Tamblyn P, Congdon C. E., 1994, ApJ, 435, 
22K 

Kroupa P, 2002, Science, 295, 82 

Leitherer C., Schaerer D., Goldader J. D, Gonzalez-Delgado 
R. M., Kune D. P., de Mello D. R, Devost D., Heckman T. P, 
1999, ApjS, 123, 3. 

Lejeune Th., CuisinierR, BuserR., 1997, A&AS, 125, 229 

Maeder A., Meynet G., 2008, RMxAC, 33, 38 

Martin-Manjon M. L., 2009, PhDT 

Martm-Manjon M. L., Molla M., Di'az A. I., Terlevich R., 2008, 

MNRAS, 385, 854 
Meynet G., Maeder A., 2005, A&AS, 429, 581 
Molla M., Garcia- Vargas M. L., 2000, A&A, 359,18 
Molla, M., Vflchez, J. M., Gavilan, M., Di'az, A. 1., 2006, MN- 
RAS, 372, 1069 

Molla M., Garcia- Vargas M. L., Bressan A., 2009, MNRAS, 398, 
451 

Moy E., Rocca-Volmerange B., Fioc M., 2001, A&A, 365, 347 
Osterbrock D. E. 1989, Astrophysics of gaseous nebulae and Ac- 
tive Galactic Nuclei. University Science Books, 1989. 
Pagel B. E. J., Edmunds M. G., Smith G., 1980, MNRAS, 193, 
219 

Panuzzo P., Bressan A., Granato G. L., Silva L., Danese L., 2003, 
A&A, 409, 99 

Panuzzo P., Silva L., Granato G. L., Bressan A., Vega, O., 2005, 
The Spectral Energy Distributions of Gas-Rich Galaxies: Con- 
fronting Models with Data, American Institute of Physics Con- 
ference Series,761.Ed.Popescu, C. C. & Tuffs, R. J. 

Pauldrach A., Hoifmann T.L., Lennon M., 2001, A&A. 375, 161 

Perez-Montero E., Diaz A. I., 2005, MNRAS, 361,1063-1076. 

Perez-Montero E., Di'az A. I., Vflchez J. M., Kehrig C, 2006, 
A&A, 449, 193 

Rauch T, 2003, A&A 403, 709 

Salpeter, E. E., 1955, ApJ, 121, 161 

Schaerer D., Vacca W. D., 1998, ApJ, 497, 618 

Silva L., Granato G. L., Bressan A., Danese L., 1998, ApJ, 509, 
103 

Smith L., Norris R., Crowther P, 2002, MNRAS, 337, 1309 
Stasinska G., 1978, A&A 66, 257 
Stasihska G., 1980, A&A, 84, 320 

Stasinska G., Collin-Souffrin S., Joly M., Alloin D., 1981, A&A, 
93, 362 

Stasihska G., 1982, A&AS, 48, 299 

Stasihska G., 1990, A&AS, 83, 501 

Stasihska G., Leitherer C, 1996, ApJSS, 107, 661-691. 

Stasihska G., 2002, RMxAC, 12, 62S. 



Stasihska G., Izotov Y, 2003, A&A, 397, 71 

Stasihska G., 2005, A&A 434,507. 

Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 

Ubeda L., Mai'z-Apellaniz J., MacKenty J. W., 2007, AJ, 133, 917 

van Zee, L., Haynes, M. P, 2006, ApJ, 636, 214V 

Vazquez G. A., Leitherer C, Schaerer D., Meynet G., Maeder A., 

2007, ApJ, 663, 995 
Vega O., Clemens M. S., Bressan A., Granato G. L., Silva L., 

Panuzzo P 2008, A&A, 484, 631 
Vflchez J. M., Pagel B. E. J., Diaz A. I., Terlevich E., Edmunds 

M. G., 1988, MNRAS, 235, 633 
Walborn N. R., Mai'z-Apellaniz J., Barba R. H., 2002, AJ, 124, 

1601 

Wotford A., 2009, MNRAS, 395, 1043 

Yin S. Y, Liang Y.C., Hammer R, Brinchmann J., Zhang B., Deng 

L. C, Plores H., 2007, A&A, 462, 535 
Zaritski D., Kennicutt R. C, Huchra J. P, 1994, ApJ, 420, 87-109 



APPENDIX A: IONIZATION PARAMETER 
RECALIBRATION. 

This appendix includes a re-calibration of the ionization parameter 
for low and high metallicity Hii regions with sulphur lines. As we 
have explained in the paper, we do not use the ionization param- 
eter as an input for our models since we derive it from the phys- 
ical properties of the ionizing clusters and the assumed geometry 
of the ionized nebula. However, since many other models use this 
magnitude to parametrize their grids or use different geometries for 
the ionized regions, we have considered interesting to obtain a re- 
calibration of the relationship between the ionization parameter and 
the sulphur lines. We have used for that the results from our models, 
as a self-consistent set of values typical for Hii regions. 

The [SlI].U6717,31/[SIII]ii9069,9532 line ratio has been 
shown to be a good ionization param eter indicator for moderate 
to high me tallicities I Diaz et alT. 
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[SII]/[SIII] to derive the ionization parameter of high metallicity re- 
gions (Z > 0.004). We have re-calibrated this relationship for each 
model. This calibration has been made for values of the ionization 
parameter logu < -2 since most Hii regions observational data do 
not show ionization parameters higher than this. Ror a more accu- 
rate calibration for higher ionization parameter values a particular 
model, according to cluster mass, should be selected a high disper- 
sion in the model results for different masses is found. 

Ror high metallicity models we find a cubic polynomial re- 
lationship between log u and log[SII]/[SIIl]. There are small dif- 
ferences between models of Z=0.008 and Z=0.004 and models of 
Z=0.02, which makes necessary to provide two different equations 
for log u in order to obtain accurate calibrations. 
For n/y = 10, we obtain; 

logu = -2.865 -1.175x +0.139x- -0.135x^ (Z=0.008, Z=0.004) 
logu = -2.705 -1.736x +0.698x^ -0.97 3x^ (Z=0.02) 
where x = log([Sll]AA671731/[SlII]AA9069,9532) 
For 11^=100, we obtain: 

logu = -3.136- O.SlOx +0.685x^ -0.784x^ (Z=0.008, Z=0.004) 
logu = -2.938 -1.051x +0.043x- -0.482x^ (Z=0.02) 
where .T = log([SlI]AA671731/[Slll]AA9069,9532) 
For low metallicity models (Z < 0.004) the calibration has 
been made using log[SII]/H„, finding a cubic polynomial relation- 
ship with log u for values of log u< -2. 
For nH=10 the fittings are: 
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Figure Al. Calibration of the ionization parameter u from tlie 
[SII]i^6717,31/[SIII]/i^9069,9532 line ratio for high metallicity models: 
Z=0.004 (black dots), Z=0.008 (red dots), Z=0.02 (green dots), and two 
particle densities, n//=10 cm"^ (upper panel) and n/f=100 cm"^(lower 
panel). The calibration has been made for values of log u< -2 and it is 
represented by orange dots. 
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Figure A2. Calibration of the ionization parameter u from the [SII]/H£, line 
ratio for the models with metallicities Z=0.0001 (black dots) and Z=0.0004 
(red dots), and both densities: n//=10cm"' (upper panel) and n//=100 cm"^ 
(lower panel). The calibration is represented by orange dots. 

logu = 18.89 +33.00.X +15.69x^ +2.33 (Z=0.0001) 
logu = -2.51 + 3.96 X +3.54 +0.69 x^ (Z=0.0004) 
where x=[SIl]AA6717,31/H„ 
For 11^=100 the fittings give the following result: 
logu = -19.95 -19.29x -7.56 x^ -1.07x^ (Z=0.0001) 
logu = -11.75 -12.98 X -6.60 X- -1.26 x^ (Z=0.0004) 
where x=[Sll]AA6717,31/H„ 



